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A huge increase of thermal conductivity k is observed at the phase transition in stoichiometric 
NaV2C>5. This anomaly decreases and gradually disappears with deviation from stoichiometry in 
Nai_ x V205 (x = 0.01, 0.02, 0.03, and 0.04). This behavior is compared with that of pure and 
Zn-doped CuGeC>3 where only modest kinks in the k(T) curves are observed at the spin-Peierls 
transition. The change of k at critical temperature T c could be partially attributed to the opening 
of an energy gap A in the magnetic excitation spectrum excluding the scattering of thermal phonons 
on spin fluctuations. However, the reason for such a strong anomaly in the k(T) may lie not only in 
the different energy scales of CuGeC>3 and NaV2C>5, but also in the different character of the phase 
transition in NaV2 0s which can have largely a structural origin, e.g. connected with the charge 
ordering. 
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Since the discovery of a spin-Peierls (SP) transition in CuGeOs [1] and NaV2 0s [2] numerous experimental and 
theoretical studies were continued to reveal the main features of this phenomenon in inorganic compounds. While 
considering that the study of CuGeOs remains very active, the main physical properties of this compound are quite 
well known by now. CuGeOs belongs to the orthorhombic space group Pbmm. The structure contains magnetic 
chains of Cu 24 " (3d 9 , S = 1/2) ions parallel to the c axis. The SP transition observed at T c ~ 14 K is marked by an 
exponential drop of the magnetic susceptibility and a simultaneous lattice dimerization along the a and c axes. Due 
to the magnetoelastic coupling the exchange constant between neighboring spins on the chain also alternates and the 
singlet ground state in the dimerized phase is separated from the band of excited triplet states by an energy gap A 
of about 25 K [3, 4]. 

In comparison with CuGeOs, only few and somewhat contradictory data are presently available on the specific 
features of the phase transition in the mixed (V 4+ /V 5+ ) or intermediate (V 4 5+ ) valence compound NaV20s. A rapid 
decrease of the magnetic susceptibility below T c ~ 35 K suggests the formation of a nonmagnetic spin-singlet ground 
state. The simultaneous formation at this temperature of a superlattice structure characteristic of a SP transition 
was confirmed by X-ray measurements [5], while the opening of an energy gap (A ~ 100 K) in the magnetic excitation 
spectrum was observed in inelastic neutron experiments [6]. These data, as well as the results of NMR [7] and ESR [8] 
studies, are strong evidence of the presence of S = 1/2 chains considered as a necessary condition for a SP transition. 
On the other hand, the two-step phase transition found in specific heat and thermal expansion measurements [9] make 
a simple SP scenario for this transition in NaV2 0s not so evident. 

The main controversy concerning the nature of the phase transition in NaV2 0s arises from the uncertainty of its 
crystal structure identification at room temperature. The X-ray diffraction pattern of single crystals fits well enough 
with both space groups P2imn and Pmmn [10-13]. In both cases the layers of edge/corner-sharing tetragonal VO5 
pyramids connected in the a — b plane are stacked along the c axis of the structure while the Na 14 " ions are situated 
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between these layers. In the former non-centrosymmetric structure, there exists alternating chains of nonmagnetic 
V 54 " ions and magnetic V 4+ ones (3d 1 , S = 1/2) aligned along the b axis (see e.g. Fig. 1 of Ref. [2]). In the latter 
structure, all the V sites are crystallographically equivalent, which results formally in a random distribution of the 
V 4+ and V 54 " ions, or in a localization of electrons on V-O-V molecular orbitals on the rungs of the ladders running 
along the b axis (see e.g. Fig. 1 of Ref. [13]). In such a situation a one-dimensional magnetic behavior does not come 
about straightforwardly but can still be justified [13-15]. 

In order to clarify a picture of the phase transition in NaV2 0s, more experimental data are necessary. One of 
the few experimental techniques not yet applied to the study of SP transitions in metal oxides is the measurement 
of thermal conductivity which can provide useful information on the interaction of elementary excitations in these 
compounds. 

In this letter we report on the first, to the best of our knowledge, such measurements, which show a surprisingly 
strong effect in NaV2 0s, in comparison with the rather weak anomaly observed in CuGeC>3. A possible interpretation 
of this effect is considered below and it is concluded that the nature of the phase transition in NaV2 0s may be more 
complicated than previously assumed. 

Single crystals of CuGeC>3 and Cuo .gsZno 02GeO3 were grown in Paris University by a floating zone method [16]. 
Single crystals of stoichiometric NaV20s were grown in Tokyo University by the self-flux method [17]. Single crys- 
tals of non-stoichiometric Nai_ :c V205 (x = 0.01, 0.02, 0.03, and 0.04) were prepared by heating a small crystal of 
stoichiometric NaV20s embedded in a large quantity of Nai_ :c V205 powder in an evacuated silica tube at 650°C for 
one week. The non-stoichiometric samples of Nai_ :c V205 contain vacancies on the Na sites. Typical dimensions of 
the crystals used in the thermal conductivity measurements were 0.2x0.8x4.0 mm 3 . The largest dimension coincided 
with the dimerization axis, i.e. the c axis in the case of CuGeOs, and the b axis in the case of NaV20s. The 
smallest dimension coincided with the a axis in the layered structure of CuGeOs and c axis in the layered structure 
of NaV2C>5. In most cases, measurements were performed in the direction of dimerization, but in CuGeOs thermal 
conductivity was measured also along the b axis. The measurements were performed by the longitudinal steady-state 
four probe method [18] with a temperature gradient equal to 2% of the current temperature. The sample holder, 
surrounded by a special shell to suppress the thermal radiation, was placed into the vacuum chamber evacuated to 
better than 10 -6 bar. The sample temperature and the temperature gradient were stabilized to better than 0.01% by 
an Oxford intelligent temperature controller with calibrated "Allen-Bradley" carbon resistors. Due to the relatively 
small dimensions of the samples the error in the determination of the absolute values of thermal conductivity was not 
better than 10%, while the precision on the thermal conductivity variations was about 1.5 %. 

The experimental results obtained in the c and b directions for pure CuGeOs and for the Zn-doped sample are 
shown in Fig. 1. A broad phonon maximum of k is observed at T = 23.3 K in the c direction, while it is not well 
pronounced in the b direction. The absolute value of k in the b direction is approximately three times lower than 
that in the c direction. The kinks on the k(T) curves are evident at the temperature of the SP transition. The SP 
transition temperature in the Zn-doped sample is reduced to T c = 10.6 K to be compared to that of the pure sample 
T c = 14.2 K. The absolute value of k at the SP transition in the doped sample is significantly smaller than that in 
the pure sample. 

The thermal conductivity of stoichiometric NaV2 0s exhibits a much more dramatic variation in the temperature 
range studied. As shown in Fig. 2, the thermal conductivity exhibits a broad maximum at T ~ 70 K and decreases 
upon approaching the SP transition. At T c a sharp upturn of k occurs so that a five times increase of thermal 
conductivity accompanied by a subsequent decrease of k upon cooling to liquid helium temperature takes place. A 
behavior similar to the one described above can be seen in non-stoichiometric samples of Nai_j;V205 (x =0.01, 0.02, 
0.03, and 0.04) with the progressively less pronounced maximumbelow T c , as shown in Fig. 3. The position of the low 
temperature maximum remains unchanged in non-stoichiometric samples but absolute values of k gradually decrease 
with x. Certain variation in the high temperature slopes and the vertical shift between the curves obtained for and 
1% and the curves for 2, 3 and 4% Na deficiency can be tentatively ascribed to a variation in the crystals morphology. 
These crystals can be easily cleaved perpendicular to the c axis, resulting in an uncontrolled extra scattering of the 
phonons on planar defects. The spin-Peierls transition temperatures in Nai_j;V205 determined from the positions 
of the minimum on the k(T) curves are shown in Fig. 4, where the solid curve is drawn as a guide for an eye. The 
suppression of T c by the Na-deficiency was observed earlier in X-ray critical scattering measurements [19]. 

Since no quantitative treatment of the transport phenomena for the spin-Peierls systems exists, we applied the 
Boltzmann equation for the analysis of the experimental data. In this approach the phonon thermal conductivity 
k ~ vlC is determined by the sound velocity v, the mean free path of phonons /, and the specific heat C which depends 
on the number of elementary excitations. The interactions of the phonon and the spin subsystems can be taken into 
account through the variation of the phonons mean free path. In the absence of these interactions, the mean free 
path of the phonons increases upon cooling and eventually saturates while their number decreases. These opposite 
tendencies result in a phonon maximum T m on the k(T) curves. The estimated phonon mean free path in CuGeOs is 
still much less than the sample dimensions; we believe that this is caused by the phonon scattering on planar defects 
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in the layered structure of this compound. The difference in the absolute values of the thermal conductivity in the 
b and the c directions in CuGeC>3 can naturally be attributed to a factor three difference in ultrasonic velocities, i.e. 
2.7- 10 3 m/sec in the b direction and 8.0T0 3 m/sec in the c direction [20]. On the other hand, the small thermal 
conductivity value of in the Cuo .gsZno 02GeO3 sample as compared with that of the pure CuGeOs sample reflects the 
increased number of scattering centers in the doped sample. 

Anomalies observed in the temperature dependence of the thermal conductivity at the spin-Peierls transition can 
be naturally connected with the change of the crystal lattice and with the formation of the energy gap in the magnetic 
excitation spectrum. The mutual positions of the temperatures of the phonon maximum T m and of the spin-Peierls 
transition T c is of crucial importance for the analysis of these anomalies. The phonon maximum in CuGeOs is reached 
at T m > T c . In this case a possible decrease of the phonon scattering on spin fluctuation at T c cannot significantly 
influence the thermal conductivity, because the phonon mean free path is already large, and only modest kinks on the 
k(T) are observed. From the behavior of NaV20s we have to conclude that the situation here is exactly opposite. If 
the maximum of k at 70 K would have the same origin, we would not be able to explain such a strong increase of k 
below T c . Therefore we have to assume that the decrease of k below 70 K in NaV20s is of a different nature, and 
most probably can be explained by the enhancement of phonon scattering on approaching the phase transition. The 
ordering occurring below T c apparently switches off this extra scattering and leads to a huge increase of the thermal 
conductivity. No variation of the sound velocity at T c can explain the observed behavior of the thermal conductivity 
in NaV2C>5. While no data for the absolute value of the sound velocity are available, a longitudinal velocity variation 
Av/v ~10 -3 was observed at T c preceded by a precursor effects observable below 70 K [21]. 

Several factors may play a role in the much stronger enhancement of the thermal conductivity below T c in NaV20s 
as compared to CuGeOs. The spin gap is much larger in NaV20s (~100 K vs ~25 K in CuGeOs). Opening of this gap 
switches off spin-phonon scattering for phonons with frequencies less that A. A four times larger spin gap in NaV20s 
implies that much more phonons will now strongly contribute to thermal conductivity below T c . This picture could, 
in principle, have allowed also to explain the behavior of the thermal conductivity of NaV20s as a function of the 
deviation from stoichiometry. Magnetic measurements show that, in contrast with the relatively small suppression 
of T c with x, the drop of magnetic susceptibility below T c is strongly reduced and for x = 0.03 - 0.04 it nearly 
vanishes. This shows that the spin gap is now, to a large extent, filled by magnetic excitations. Correspondingly, the 
phonon-spin fluctuation scattering is being restored, which suppresses the increase of k below T c in non-stoichiometric 
samples. 

However, the comparison with CuGeOs shows that pure spin-phonon scattering is hardly sufficient to explain huge 
increase of k observed. Magnetoelastic coupling in NaV20s is not strong enough, which follows e.g. from the very 
weak dependence of T c on magnetic field [23]. We believe that the huge increase of the thermal conductivity in 
NaV205 can be related to the specific features of the phase transition in this compound. If the crystal structure 
above T c is indeed that obtained in recent studies [11-13], all V ions are on the average equivalent, i.e. the V 4 5+ 
ions. Such a state can be visualized as the (dynamic) random mixture of V 44 " and V 54 ". If slow enough, these charge 
fluctuations can strongly enhance phonon scattering and reduce thermal conductivity similar to a situation in glasses. 
The phase transition may then be accompanied, or even caused, by the charge ordering [22] at which this extra 
scattering mechanism is switched off. The opening of a spin gap may then be not so much a driving force but a 
consequence of this ordering. The transition from the disordered glass-like phase to an ordered one could explain the 
strong enhancement of the thermal conductivity below T c . This picture can also explain the observation [23] that the 
total entropy of the transition in NaV20s is larger than the pure spin one. 

Although our measurements do not allow to draw a final conclusion about the microscopic picture of the changes 
occurring at the phase transition in NaV2 0s, the qualitative difference observed in the behavior of CuGeOs and 
NaV2 05 can be taken as an evidence that the nature of the transition in NaV2 0s is more complicated than in 
CuGeOs and most probably directly involves lattice and charge degrees of freedom, e.g. in a form of charge ordering. 
Further experiments, especially neutron scattering, should shed light on the microscopic picture of the spin-Peierls 
transition in NaV20s and should clarify the nature of the strong anomaly in the thermal conductivity of NaV20s 
reported in this letter. 
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FIG. 1. Temperature dependences of thermal conductivity in pure and Zn-doped CuGeOs. 

FIG. 2. Temperature dependences of thermal conductivity k and of magnetic susceptibility x m stoichiometric 
NaV 2 5 . 

FIG. 3. Temperature dependences of thermal conductivity in Nai_ :c V205 (x = 0, 0.01, 0.02, 0.03, and 0.04). 
FIG. 4. Critical temperatures of the spin-Peierls transition in Nai_ :c V205 as obtained from the thermal conductivity 
data. The solid curve is drawn as a guide for an eye. 
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